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We present new millimeter and radio observations of nine z~6 quasars dis- 
covered in deep optical and near-infrared surveys. We observed the 250 GHz 
r|j ■ continuum in eight of the nine objects and detected three of them. New 1.4 GHz 

(-H ' radio continuum data have been obtained for four sources, and one has been de- 

Qh' tected. We searched for molecular CO (6-5) line emission in the three 250 GHz 

O ■ detections and detected two of them. Combining with previous millimeter and 
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radio observations, we study the FIR and radio emission and quasar-host galaxy 
evolution with a sample of 18 z~6 quasars that are faint at UV and optical 
wavelengths (rest-frame 1450 A magnitudes of mi45o > 20.2). The average FIR- 



^ ' to-AGN UV luminosity ratio of this faint quasar sample is about two times higher 



than that of the bright quasars at z~6 (mi45o < 20.2). A fit to the average FIR 
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O ■ and AGN bolometric luminosities of both the UV/optically faint and bright z~6 
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quasars, and the average luminosities of samples of submillimeter/millimeter- 
observed quasars at z~2 to 5, yields a relationship of LpiR ~ Lboi'^'^^- Five of the 
18 faint z~6 quasars have been detected at 250 GHz. These 250 GHz detections, 
as well as most of the millimeter-detected optically bright z~6 quasars, follow a 
shallower trend of Lpm ~ Lboi'^'^^ defined by the starburst-AGN systems in local 
and high-z universe. The millimeter continuum detections in the five objects and 
molecular CO detections in three of them reveal a few xlO® Mq of FIR-emitting 
warm dust and 10^° Mq of molecular gas in the quasar host galaxies. All these 
results argue for massive star formation in the quasar host galaxies, with esti- 
mated star formation rates of a few hundred M0yr~^. Additionally, the higher 
FIR-to-AGN luminosity ratio found in these 250 GHz-detected faint quasars also 
suggests a higher ratio between star formation rate and supermassive black hole 
accretion rate than the UV/optically most luminous quasars at z~6. 

Subject headings: galaxies: quasars — galaxies: high-redshift — galaxies: star- 
burst — molecular data — galaxies: active 



1. Introduction 

Understanding the formation and evolution of the first galaxies and the supermassive 
black holes (SMBHs) they host is a key goal of both observational and theoretical astronomy. 
The first sample of more than twenty bright quasars at z~6 were selected from the imaging 
data of the Sloan Digital Sky Survey (hereafter SDSS main survey, e.g.. Fan et al. 2000, 
2006). These objects are very bright in quasar UV and optical emission (SDSS magnitudes 
of 18.74 < zab < 20.42) and have SMBHs with masses of Mbh ~ a few 10^ Mq (Jiang 
et al. 2007; Kurk et al. 2007). The dust and molecular gas properties of this bright 
quasar sample at z~6 have been studied at millimeter wavelengths, using the Max Planck 
Millimeter Bolometer Array (MAMBO), the IRAM Plateau de Bure Interferometer (PdBI), 
and the Very Large Array (VIA, Bertoldi et al. 2003a; 2003b; Petric et al. 2003; Walter et 
al. 2003, 2004, 2009; Carilh et al. 2004, 2007; Riechers et al. 2009; Wang et al. 2007, 2008b, 
2010). About 30% of them have strong continuum emission at 250 GHz from 40-to-60 K 
warm dust (Bertoldi et al. 2003a; Petric et al. 2003; Wang et al. 2007, 2008b). The FIR 
luminosities of the MAMBO detections are a few xIO^^Lq to 10^^ Lq, and the estimated 
dust masses are a few xlO^M© (Beelen et al. 2006; Bertoldi et al. 2003a; Wang et al. 
2008a). Most of the 250 GHz-detected sources have also been detected in molecular CO 
line emission with PdBI, yielding molecular gas masses of (0.7 — 3) x 10^° M© (Walter et al. 
2003; Carilh et al. 2007; Riechers et al. 2009; Wang et al. 2010, 2011). The FIR-CO and 
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FIR-radio luminosity ratios of these objects are consistent with the trend defined by dusty 
starburst systems at low and high redshifts (Riechers et al. 2006; Beelen et al. 2006; Wang 
et al. 2008b). The results provide direct evidence for massive star formation at a rate of 
~ 10^ Mq yr~^, indicating that bulges are being formed at the same time that the SMBH is 
rapidly accreting in at least 1/3 of the optically bright quasar systems at z~6 (Wang et al. 
2008b). Molecular CO (3-2) and (7-6) line emission was resolved in one of the FIR and CO 
luminous quasars, SDSS J114816. 64+525150. 3 at z=6.42, revealing a source size of ~ 5kpc 
(Walter et al. 2004, 2009). The dynamical mass estimated with the size and line width of 
the CO emission suggest a SMBH- spheroidal bulge mass ratio an order of magnitude higher 
than the typical present-day value (Walter et al. 2004). 

Z~6 quasars that are fainter at UV and optical wavelengths have been discovered from 
deep optical-near infrared surveys (the SDSS southern deep imaging survey, 20.54 < zab < 
22.16, Jiang et al. 2008, 2009, the Canada- France High redshift Quasar Survey [CFHQS], 
20.26 < ZAB < 24.40, Willott et al. 2007, 2009a, 2009b). Near-infrared spectroscopic 
observations of some of these faint quasars reveal that they are likely to be accreting at the 
Eddington rate, and have less massive black holes (i.e., 10® to 10^ M©) than the SDSS main 
survey z~6 quasar sample (Kurk et al. 2009; Willott et al. 2010b). The steepness of the 
quasar luminosity function at high redshifts is such that lower-luminosity/mass objects are 
more common. Thus, these faint z~6 quasars provide an important sample to understand 
SMBH-host galaxy formation at the earliest epoch. 

The millimeter dust continuum emission from the faint z~6 quasars was first studied 
with a small sample of four objects discovered by CFHQS (Willott et al. 2007). One of 
them was detected at 250 GHz with a fiux density of ~ 1 mJy, indicating a FIR continuum 
luminosity from dust of > 10^^ Lq. The mean 250 GHz fiux density of the four quasars when 
stacked is ~ 0.6 mJy, which is about a factor of two smaller than the average emission of 
the z~6 luminous quasar sample (Willott et al. 2007; Wang et al. 2008b). We subsequently 
detected 250 GHz dust continuum and molecular CO (6-5) line emission in another two z~6 
quasars with zab ~ 20.7 (SDSS J205406.42-000514.8 and NDWFS J142516.30+325409.0, 
Jiang et al. 2008; Cool et al. 2006; Wang et al. 2008b, 2010) and CO (2-1) line emission 
from one CFHQS quasar at z=6.2 (CFHQS J142952+5447170, Willott et al. 2010a; Wang et 
al. 2011). These observations revealed a few 10® M© of 40-to-60 K warm dust and 10^° Mq 
of molecular gas in the quasar host galaxies, i.e., masses comparable to those found with the 
bright z~6 quasars. 



-'^We exclude this object from our analysis throughout this paper as there is no published millimeter 
continuum observation yet. 
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In this paper, we report new millimeter and radio observations of nine z~6 quasars. We 
then study the FIR and radio emission with a sample of 18 z~6 quasars that are faint at UV 
and optical wavelengths, and investigate the dust and molecular gas emission properties and 
star forming activity in the host galaxies of the millimeter-detected objects. We summarize 
the current quasar sample at z~6 and describe the new millimeter and radio observations in 
Section 2, present the results in Section 3, analyse the continuum emission properties of the 
faint z~6 quasar sample in Section 4, discuss the star formation and quasar-host evolution 
in Section 5, and conclude in Section 6. A A-CDM cosmology with Hq = 71km s^^ Mpc^^, 
Qm = 0.27 and Q/^ = 0.73 is adopted throughout this paper (Spergel et al. 2007). 



2. Sample and observations 

2.1. Sample 

There are thirty-three z~6 quasars that have published millimeter and radio observa- 
tions from the literature (Bertoldi et al. 2003a; Petric et al. 2003; Carilli et al. 2004; Wang 
et al. 2007; 2008b; Willott et al. 2007). In this paper, we present new MAMBO observations 
of the 250 GHz dust continuum emission from eight z~6 quasars (Table 1), and molecular 
CO observations of three of them. Seven of the eight sources were discovered from the SDSS 
southern survey with optical magnitudes of 20.87 < zab < 22.16 (Jiang et al. 2008, 2009), 
and the other source, ULAS J131911.29+095051.4 (hereafter J1319+0950), was selected from 
the UKIRT Infrared Deep Sky Survey with zab = 19.99 (UKIDSS, Mortlock et al. 2008). 
Seven of the eight objects (all but SDSS J020332.39+001229.3, hereafter J0203+0012) have 
never been observed at millimeter wavelengths before. We also report new radio 1.4 GHz 
continuum observations for three of eight new MAMBO-observed objects and another z~6 
quasar from CFHQS (Willott et al. 2007). 

The new observations, together with the previous data, now provide us a sample of 
forty 250 GHz-observed z~6 quasars with optical magnitudes over a wide range, i.e., 18.74 < 
Zab < 22.16, or quasar rest-frame 1450 A AB magnitudes of 18.81 < mi45o < 22.28. Wang 
et al. (2008b) studied the dust emission from quasars at z~6, mainly based on the most 
luminous sample from the SDSS main survey. In this paper, we investigate the FIR and 
radio emission properties of the fainter quasars at z~6, i.e., a subsample of 18 objects with 
mi45o > 20.2 (Table 1 and 2), and compare them to the bright z~6 quasars (mi45o < 20.2) 
discussed in previous papers (Bertoldi et al. 2003a; Petric et al. 2003; Wang et al. 2008b). 
We will discuss the optical measurements and quasar UV and bolometric luminosities based 
on mi45o, which is derived from the SDSS z-band magnitudes and deep optical spectra (e.g.. 
Fan et al. 2000, 2001, 2003, 2004, 2006; Jiang et al. 2008, 2009; Willott et al. 2007), or 
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from infrared photometry in the cases in which significant dust reddening is presented in the 
optical spectra (e.g., McGreer et al. 2006). 



2.2. Observations 

Observations of the 250 GHz dust continuum from the eight z~6 quasars were carried out 
using the MAMBO-II 117-element array on the IRAM 30-m telescope (Kreysa et al. 1998) in 
the winters of 2008-2009 and 2009-2010. We adopted the standard on-off photometry mode 
with a chopping rate of 2 Hz by 32" in azimuth (see also Wang et al. 2007, 2008). The on- 
source observing time was about 1-3 hours for each object, to achieve la rms sensitivities of 
~ 0.5 to 0.7 mJy. We used the MOPSIC pipehne (Zylka 1998) to reduce the data, and three 
sources, J1319+0950, SDSS J012958.51-003539.7 (hereafter J0129-0035), and J0203+0012, 
were detected at > 4a. Wang et al. (2008b) observed J0203+0012 with MAMBO and did 
not detect the source with a la rms of 0.7 mJy. Our new observations, together with the old 
data, improve the rms to 0.46 mJy for this object, and we detect the 250 GHz continuum at 
4a. We have also obtained a 2.4o" signal from another source, SDSS J023930. 24— 004505.4 
(hereafter J0239-0045). 

We also observed four z~6 quasars at 1.4 GHz in Winter 2008-2009, using the VLA 
in A configuration with a synthesized beam size of FWHM ~ 1.4". Flux scales were mea- 
sured using the standard VLA calibrators of 3C286 and 3C48, and phase calibrations were 
performed every 30 min on nearby phase calibrators (1254-116, 1507-168, 2136-1-006, and 
2323-032). The typical rms noise level was ~ 17/iJybeam"^ in three hours of observing 
time for each source. We reduced the data with the standard VLA wide field data reduction 
software AIPS, and one source, J1319-I-0950, was detected at ~ 3.8o". 

We searched for molecular CO (6-5) line emission from the three MAMBO detections 
using the 3 mm receiver on the PdBI. The source J1319-I-0950 was observed in Summer 2009 
in C and D configurations (i.e. a spatial resolution of FWHM ~ 3.5") with the narrow-band 
correlator, which provides a bandwidth of 1 GHz in dual polarization. The observations of 
J0129— 0035 and J0203-I-0012 were carried out using the new wide band correlator WideX 
with a bandwidth of 3.6 GHz in dual polarization, which covers well the Lya-redshift un- 
certainty of Az ~ ±0.03 (Jiang et al. 2009). The observations were performed in Summer 
2010 in D configuration with a synthesized beam size of FWHM ~ 5". Phase calibration was 
performed every 20 min with observations of phase calibrators, 1307+121 and 0106-1-013, 
and we observed MWC 349 as the fiux calibrator. We reduced the data with the IRAM 
GILDAS software package (Guilloteau & Lucas 2000). The CO (6-5) line emission is de- 
tected from J1319+0950 at 4. 80" in 10.3 hours of on-source time with an rms sensitivity of 
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0.50mJybeam~^ per 70kms~^ binned channel, and from J0129— 0035 at 5.3a in 8.8 hours 
with an rms of 0.55 mJy beam" ^ per lOOkms"^ channel. We also detect the continuum at 
> 3(T at the CO (6-5) line frequency in both objects. The other source, J0203-I-0012, was un- 
detected in both CO and 3 mm continuum in 13.9 hours of on-source time, with a sensitivity 
of 0.44mJybeam~^ per 100 km s"^ channel. 



3. Results 

3.1. Notes for individual objects 

We summarize the new millimeter and radio observations of the nine z~6 quasars in 
Table 1. Detailed measurements of the three new 250 GHz dust continuum/molecular CO 
(6-5) detections are listed below: 

J1319-|-0950 The source was discovered in the UKIRT Infrared Deep Sky Survey 
(UKIDSS) with mi45o = 19.65, which lies in the typical magnitude range of the opti- 
cally bright z~6 quasars selected from the SDSS main survey (Mortlock et al. 2008). 
We detect the 250 GHz dust continuum emission from this object with a flux density of 
S250GHZ = 4.20 ± 0.65 mJy, which is the third-strongest 250 GHz detection among the known 
z~6 quasars. The 1.4 GHz radio continuum is also detected from our VLA observation, 
with Si.4GHz = 64 ± 17/iJy (left panel of Figure 1). The radio emission appears unresolved 
with a peak at 13^^19^11. 30% 09°50'51.8", i.e. within 0.4" of the optical quasar position. 
We also see a Aa peak (Si.4ghz = 68 ± 17/iJy) at 13^19"'11.65% 09°50'53.2" on the ra- 
dio map, which is 5.6" away from the quasar. For comparison, only 0.02 detections with 
Si.4GHz > 60/iJy are expected from a random circular sky area with a radius of 6", ac- 
cording to the 1.4 GHz radio source counts from previous deep VLA observations (i.e., 
N(> Si.4GHz) = (0.40±0.04)(Si.4GHz/75/iJy)(-i-43±o.i3)^ Fomalont et al. 2006). No optical 
counterpart has been found at the radio source position in the NASA/IPAC Extragalactic 
Database or the SDSS images. 

We observed the CO (6-5) line emission from this object with the PdBI. The line is 
detected at 4.8o" on the velocity-integrated map with a line flux of 0.43 ± 0.09 Jykms^^, 
0.9" away from (i.e., consistent with) the quasar position (the upper panel of Figure 2). A 
Gaussian fit to the line spectrum yields a peak line flux density of 0.80±0.15 mJy beam~^ and 
a line width of FWHM = 537 ± 123kms-\ centered at the redshift of z = 6.1321 ± 0.0012. 
This is slightly higher than, but consistent with the Mg II redshift of z = 6.127 ± 0.007, 
considering the measurement errors (±0.004, Mortlock et al. 2008) and the typical velocity 
shift and scatters 97 ± 269 km s^^ (—0.002 ± 0.006 in redshift) between the broad Mg II 
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line emission and the quasar systemic redshift (Richards et al. 2002). We also detect 
the continuum under the CO (6-5) hne emission with a peak value of SgycHz = 0.31 ± 
0.08mJybeam~^ on the continuum map averaged over the line-free channels (right panel of 
Figure 1). We adopt this as the continuum flux density as the source is unresolved with the 
5" beam. We checked the position of the 1.4 GHz radio source 5.6" away from the quasar in 
the PdBI CO/97 GHz continuum images. There is no clear (> 3cr) detection at this position 
and the measurement in the 97 GHz continuum image is 0.19 ± 0.08 mJy (right panel of 
Figure 1). 

J0129— 0035 This is the faintest of the SDSS z~6 quasars, with mi45o = 22.16 (Jiang 
et al. 2009). We detect 250 GHz dust continuum from the quasar host galaxy, with a flux 
density of S250GHZ = 2.37 ± 0.49 mJy, comparable to the flux densities of the optically bright 
z~6 quasars with 250 GHz detections. 

We detected the CO (6-5) line in this source with a flux of 0.37 ± 0.07 Jykms^^, i.e., 
at 5.3cr, integrated over a velocity range of 380 kms~^ (the middle panel of Figure 2). The 
line width fltted with a single Gaussian is 283 ± 87kms^^, with a line peak of 1.2±0.3 
mJy. The host galaxy redshift measured with the CO line is z = 5.7794 ± 0.0008, which 
is in good agreement with the optical redshift of ZLya = 5.78 ± 0.03 measured with Lya 
(Jiang et al. 2009). The underlying dust continuum is also detected with a flux density of 
Sio2GHz = 0.14 ± 0.04 mJy. 

J0203-|-0012 This object is a broad absorption line (BAL) quasar discovered by both 
the SDSS southern survey (Jiang et al. 2008) and the UKIDSS survey (Venemans et al. 
2007; Mortlock et al. 2008) with mi45o = 20.97. The redshift measured with the quasar 
NV line emission is z=5.72 (Mortlock et al. 2008). Wang et al. (2008b) detected the 1.4 
GHz radio continuum with Si.4ghz = 195 ± 22/iJy (Wang et al. 2008b). Our MAMBO 
observation yields a 250 GHz dust continuum flux density of S250GHZ = 1-85 ± 0.46 mJy. We 
have searched for the CO (6-5) line emission from this object over a wide redshift range, from 
z=5.60 to 5.84, but did not detect it in 14 hours of on-source time. Assuming a line width 
of 600 kms~^, the 3cr upper limit of the CO flux is estimated to be 0.32 Jykms"^ (the lower 
panel of Figure 2p. The dust continuum at the CO (6-5) line frequency is also undetected 
with a continuum sensitivity of 0.04mJybeam~^ and a 3cr upper limit of 0.12 mJy. 



^We estimate the 3a upper limit for the CO line flux as 3crchannci(AViinGAVchannci)^^^, where AVunc is 
the assumed line width, and iTchannci and AVchannci are the rms per channel and channel width listed in 



Section 2.2 (Seaquist et al. 1995). 



3.2. FIR and radio luminosities and dust mass 

We calculate the FIR luminosities and upper limits for the eight z~6 quasars observed 
by MAMBO. We model the FIR emission with an optically thin graybody and normalize the 
model to the MAMBO-250 GHz and available PdBI-3mm dust continuum flux densities/So" 
upper limits (see also Wang et al. 2007, 2008). A dust temperature of Td = 47 K and 
emissivity index of /3 = 1.6 are adopted in the model; these are typical values for the 
FIR and CO luminous quasars at lower redshifts (Beelen et al. 2006). We then integrate 
the model FIR SED from 42.5 /im to 122.5 /im to calculate the total FIR luminosities/upper 
limits. The FIR luminosity of the strongest MAMBO detection, J1319+0950, is 1.0 x 10^^ L^ 



05 



making it one of the most FIR-luminous objects among the z~6 quasar sample. The other 
two MAMBO-detected objects are from the SDSS southern survey sample and have FIR 
luminosities of 4 ~ 5 x 10^^ L0. We also estimate the dust masses for the three MAMBO 
detections, using Mdust = LFiR/47r J K^Q^dv, where Bj^ is the Planck function and k^j = 
^^{v/vqY is the dust absorption coefficient with kq = 18.75 cm^g^-*^ at 125 ;um (Hildebrand 
1983). The estimated dust masses of the three sources are (2.5 — 5.7) x 10® Mq. We finally 
calculate the 1.4 GHz radio continuum luminosity densities and upper limits for the four 
VLA observed objects. We assume a power-law continuum of i,y ~ ^-0.75 (^Qondon 1992) 
to correct the observed frequency to the quasar rest frame. The results are summarized 
in Table 3; we also include the 1.4 GHz radio luminosity of J0203+0012 from Wang et al. 
(2008b). 



3.3. CO luminosities and molecular gas masses 

We calculate the CO (6-5) line luminosities and upper limits for the three CO-observed 
z~6 quasars, using Lco(6-5) — 3-25 x 10^IAvco(6-5)'^obs~^DL^(l + z)~^ (Kkms""*^ pc^), where 
IAvco(6-5) is the line fiux or 3a upper limit in Jykms"^ presented in Table 1, z/obs is the 
observed frequency of the CO (6-5) line in GHz, Dl is the luminosity distance in Mpc, and z is 
the quasar redshift (Solomon & Vanden Bout 2005). We then assume a CO excitation ladder 
similar to that of J1148-I-5251 (Riechers et al. 2009), for which L'^^Q/g^-j/L'^^Q^Q-j ^ 0.78, 
allowing us to calculate the CO (1-0) luminosities. The FIR-to-CO(l-O) luminosity ratios 
of J1319+0950 and J0129-0035 are 520 and 340 Lg/Kkms"^ pc^, respectively, which are 
similar to the values found for other z~6 CO-detected quasars (Wang et al. 2010). We then 
calculate the molecular gas masses (Mgas = M[H2 + He]), adopting a CO luminosity-to-gas 
mass conversion factor of a = O.8M0 (Kkms"^pc^)~^ seen in low-redshift ultra-luminous 
infrared galaxies (ULIRGs, Solomon et al. 1997; Downes & Solomon 1998). The molecular 
gas masses for the two new CO detections are all of order 10^° Mq (Table 3). 
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4. Analysis 

The new observations increase the sample of MAMBO-observed z~6 quasars to forty 
objects (Bertoldi et al. 2003a; Petric et al. 2003; Willott et al. 2007; Wang et al. 2007, 
2008b). In this paper, we study their FIR and radio emission properties, focusing on the 18 
sources that are fainter at UV and optical wavelengths (20.20 < mi45o < 22.28, hereafter 
faint z~6 quasars) and compare them with the first sample of millimeter-studied z~6 quasars 
from the SDSS main survey (hereafter bright z~6 quasars, Bertoldi et al. 2003a; Petric et al. 
2003; Wang et al. 2007, 2008b). The MAMBO observations of these faint z~6 quasars have 
typical sensitivities (0.4 to 0.6 mJy) comparable to or slightly lower than that achieved in 
previous MAMBO observations of the SDSS main survey z~6 quasars (0.5 — 1.1 mJy, Petric 
et al. 2003; Wang et al. 2007, 2008) and the optically bright quasar samples at redshifts 
2 and 4 (Omont et al. 2001,2003; Carilh et al. 2001). Five of the 18 sources are detected 
at 250 GHz at > 3ct (Wang et al. 2007, 2008b; Willott et al. 2007), yielding a detection 
rate of 28±12%, or 25 ± 14% if we consider the subsample of twelve objects from the SDSS 
southern survey selected with the same color selection criteria (three detections. Fan et al. 
2004; Jiang et al. 2008; 2009; Wang et al. 2008b, this work). These are slightly lower than 
the MAMBO detection rate of >30% found with the optically bright quasar samples from 
z=2 to 6 (Omont et al. 2001, 2003; Carilh et al. 2001; Wang et al. 2008b), but are consistent 
within the errors. 

We searched for molecular CO (6-5) line emission in four of the five MAMBO-detected 
faint z~6 quasars, and detected three (Wang et al. 2010), indicating that highly excited 
molecular gas is ubiquitous in these millimeter (FIR) luminous quasar-galaxy systems at 
the earliest epoch. Fourteen of the 18 sources have radio 1.4 GHz continuum observations 
with typical la rms of ~ 20 /xJy, and only two of them, J0203+0012 and J1427+3313j have 
radio detections (McGreer et al. 2006; Wang et al. 2007, 2008b, this work). This gives 
a much lower radio detection rate for the faint z~6 quasar sample than that of the bright 
objects; nine out of 22 z~6 quasars with mi45o < 20.2 were detected at 1.4 GHz with similar 
observational sensitivities (~ 20/iJy, Petric et al. 2003; Carilli et al. 2004; Wang et al. 
2007, 2008b). The 1.4 GHz continuum flux density is 1.7 mJy for J1427+3312 (McGreer et 
al. 2006) and 195 /iJy for J0203+0012 (Wang et al. 2008b), indicating radio-loud nuclei in 
the two objectqj (Momjian et al. 2008). Deep observations are required to better constrain 



^Note that J1427+3312 was selected from radio data (McGreer et al. 2006). 

^The radio-loud quasars are defined as quasars with radio-to-optical flux density ratios of 
fj^5GHz/f ° > 10 (Kellermann et al. 1989) where fluxes are measured in the rest frame. We adopt 

this definition and calculate R assuming a power-law UV-to-optical continuum of i^, ^ v^'^-^ (Fan et al. 
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the radio emission from the radio-quiet UV/optical faint z~6 quasars. 



4.1. The average FIR and radio emission 

We calculate the average 250 GHz and 1.4 GHz continuum for the 18 faint z~6 quasars. 
The mean flux densities are derived for i) all sources, ii) MAMBO detections, and iii) 
MAMBO non-detections, using the MAMBO and VLA measurements weighted by inverse 
variance (See also Wang et al. 2008b). We also present the average values for the 22 bright 
quasars (mi45o < 20.2) for comparison. The average emission of this bright quasar sample is 
consistent with the results presented in Wang et al. (2008b). We list all the results in Table 
4. 

The mean 250 GHz flux density is 0.85 ± 0.12 mJy for all 18 faint quasars, and 1.85 ± 
0.20 mJy for the five MAMBO-detected sources, which are about 50% and 70% of the values 
found for the bright quasars, respectively. The results are about the same when the two 
radio-loud quasars are excluded. The stacking measurement for the 13 MAMBO-undetected 
sources in the faint sample is 0.33 ± 0.15 mJy, which places an upper limit of < 0.78 mJy 
(stacking average plus the Scr rms) for the average millimeter dust continuum emission. 

We convert the average 250 GHz flux densities to 42.5 122. 5 //.m FIR luminosities, 

adopting the assumptions described in Section 3.2, and calculate the average FIR-to-1450A 
luminosity ratios ((Lfir) / (L1450)) for each sample (see Table 4). The average luminosity 
ratios for all the sources and the MAMBO detections in the faint-quasar sample are two 
and three times higher than the corresponding values found with bright quasars. This is 
consistent with the non-linear FIR-to-AGN luminosity relationships of Lfir ~ L^^p^, where 
a < 1 found with the AGN-starburst at low and high redshifts (Beelen et al. 2004; Hao et 
al. 2005, 2008; Wang et al. 2008b; Serjeant & Hatziminaoglou 2009; Lutz et al. 2010), as 
we will discuss further in Section 4.2. However, the different average luminosity ratios found 
with the MAMBO detections in the bright and faint quasar samples may be partly due to 
the selection effect that the detection limit of the MAMBO observations are about the same 
for the two samples, while the AGN luminosities are different by design. 

The average 1.4 GHz continuum fiux density, derived from 14 VLA-observed sources 
in the faint quasar sample, is 19 ± 6 /iJy, and 2 ± 6 /iJy when the two radio-loud sources 
are excluded. This suggests an upper limit of 20 /iJy for the radio-quiet and UV/optically 
faint z~6 quasars. We estimate the average rest-frame 1.4 GHz luminosities/upper limits 



2000) and radio continuum of f,y ^ i^ °-^^ (Condon 1992). 
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for the two samples and list the results in Table 4. We plot Lpm vs. Li^^i.4ghz of all the 
MAMBO-detected z~6 quasars together with the average luminosities of the MAMBO- 
detections and non-detections in Figure 3, and compare them to sample of IRAS-selected 
local star forming galaxies (Yun et al. 2001). The local star forming galaxies with FIR 
and radio emission powered by the active star formation have an average luminosity ratio 
of q = log (Lfir/3.75 x 10^^ l^) _ i^g (Li.4ghz/L0 ■ Hz"^) = 2.34. The range 1.6 < q < 3.0 
contains 98% of the IRAS galaxies (Yun et al. 2001). The average luminosities/upper limits 
of the MAMBO-detected z~6 quasars, though lower than the mean value of q = 2.34, are 
all within this range defined by local star forming galaxies. This argues that, in addition to 
the AGN power, star formation also contributes significantly to the FIR-to-radio emission 
of the MAMBO-detected z~6 quasars. The average ratio of the MAMBO-undetected bright 
z~6 quasars is offset from this range. 



4.2. The FIR-to-AGN luminosity relationship 

The correlations between FIR luminosities and AGN UV, optical, and bolometric lumi- 
nosities have been studied widely with samples of infrared and optically selected quasars to 
understand the origin of quasar FIR emission. Study of samples of the Palomar-Green (PG) 
quasars and Seyfert galaxies shows a relationship of LpiR ~ L^q^°'^ for typical optically- 
selected AGN/quasars in the local universe (Hao et al. 2005; Netzer et al. 2007; Lutz et 
al. 2010). A shallower relationship of LpiR ~ ^%gn ^^.s found for the local IR-luminous 
quasars hosted by starburst ULIRGs (Hao et al. 2005) and high-redshift FIR-and-CO lu- 
minous quasars (Beelen et al. 2004; Hao et al. 2008; Wang et al. 2007; Lutz et al. 2010), 
which suggests excess FIR emission powered by extreme starburst in these systems. Wang 
et al. (2008b) reported that the millimeter-detected bright z~6 quasars follow the shallow 
luminosity trend defined by the local IR quasars. 

In this work, we investigate the FIR-to-AGN luminosity correlation of the full sample 
z~6 quasars, and compare with the samples of MAMBO-250 GHz or SCUBA-350 GHz ob- 
served optically bright quasars at z~2 to 5 (Omont et al. 2001, 2003; Carilli et al. 2001; 
Priddey et al. 2003a), and local IR quasars (Zheng et al. 2002; Hao et al. 2005). We convert 
the UV or optical luminosities of all objects to AGN bolometric luminosities assuming a UV- 
to-optical power-law continuum of ii, ~ z/""'^, and an optical-to-AGN bolometric luminosity 
conversion factor of Lboi = 10.4 z/L ° (Richards et al. 2006). In figure 4, We plot all 

the z~6 quasars, as well as the local IR quasar sample and the average luminosities of the 
submillimeter/millimeter [(sub)mm]-observed z~2 to 5 quasar samples (Omont et al. 2003; 
Hao et al. 2008). The average FIR luminosities of the z~2 to 5 quasars are derived with the 
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mean flux densities at 250 GHz or 350 GHz (Omont et al. 2001, 2003; Carilli et al. 2001; 
Priddey et al. 2003a), adopting the same assumptions of dust temperature and emissivity 
index used for the z~6 quasars (See Section 3.2). We fit a fine to the average luminosities 
of all these (sub)mm-observed quasar samples at z~2 to 6 using the Ordinary Least Square 
method (Isobe et al. 1990). This yields an relationship of: 

log (y^) = (0-62 ± 0.09)log (^\ + (3.90 ± 1.29) (1) 

between the average FIR and AGN emission in these high-redshift optically-selected quasars, 
and the non-linear slope of 0.62 might suggest contributions from both host galaxy star 
formation and AGN power to the FIR emission (Beelen et al. 2004; Netzer et al. 2007; Hao 
et al. 2005, 2008; Wang et al. 2008b; Serjeant & Hatziminaoglou 2009; Lutz et al. 2010) 

The MAMBO-250 GHz observations at sub-mJy sensitivity have detected the most FIR 
luminous objects among the z~6 quasars (Bertoldi et al. 2003a, Petric et al. 2003; Wang et 
al. 2007, 2008, this work). Most of these objects follow the shallower luminosity correlation 
trend of the local IR quasars and the (sub)mm-detected quasars at z~2 to 5. In particular, 
the five MAMBO-detected z~6 quasars with m^so > 20.2 (Willott et al. 2007; Wang et al. 
2008b; this work) overlap with the high luminosity end of the IR quasars. This may suggest 
a starburst-dominant FIR emission in these most FIR luminous quasars at z~6, similar to 
that found in the local IR quasars. A fit to all the (sub) mm detections in the high-z quasar 
sample and the local IR quasars gives a relationship of 

log ( Y^) = (0-45 ± 0.03)log (j^\ + (6.62 ± 0.39) (2) 



5. Discussion: star formation in the MAMBO-detected faint z~6 quasars 

Strong FIR dust continuum emission has been detected in the host galaxies of >l/4 
of the faint z~6 quasars (mi45o > 20.2, Willott et al. 2007; Wang et al. 2008b, 2010, 
this work). The dust continuum and molecular CO detections reveal a few xIO^Mq of 
FIR-emitting warm dust and 10^'^ Mq of molecular gas in the host galaxies of these faint 
quasars, which are comparable to the typical values found with the millimeter-detected 
optically bright SDSS z~6 quasars. The dust formation process associated with the evolution 
of massive stars on time scales < IGyr has been discussed widely since the discovery of 
4 X 10^ M0 of dust in one of the most luminous/massive z~6 quasars, J1148+5251 (Beelen 
et al. 2006; Morgan & Edmunds 2003; Maiolino et al. 2004; Dwek et al. 2007). The 
detections of comparble amounts of dust in these faint z~6 quasars may suggest similar fast 
dust formation mechanism in their host galaxies. Their FIR luminosities estimated from the 
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millimeter measurements are 3 to 5 x 10^^ Lq, falling on the FIR-to-AGN luminosity trend 
defined by the local IR quasars and the FIR luminous quasars at z~2 to 5 (Section 4.2). 
The lower limit of their average FIR-to-radio emission ratio of q > 1.6 derived with the 250 
GHz and 1.4 GHz measurements is within the typical range of 1.6 < q < 3.0 in star forming 
galaxies (Section 4.1). 

All these results suggest significant star formation in the millimeter-detected optically 
faint quasars at the highest redshift. Conservatively assuming that 50% of the FIR emission 
is powered by host galaxy star formation, the average star formation rate (SFR) of the five 
millimeter-detected mi45o > 20.2 quasars is about 560 M0yr~^ (Kennicutt 1998). Three 
of them have been detected in molecular CO (6-5) line emission (Wang et al. 2010; this 
work) and the ratios between FIR and CO luminosities are 220 to 370 LQ/Kkms~^ pc^, 
which is within the typical range found with the submillimeter galaxies and CO-detected 
quasars at lower redshifts (Solomon & Vanden Bout 2005; Greve et al. 2005). This result 
suggests a high star formation efficiency in these CO-detected faint z~6 quasars, similar to 
that of the extreme starburst systems (Solomon & Vanden Bout 2005). Additionally, if the 
MAMBO-detected faint z~6 quasars are indeed less massive SMBH systems accreting at 
their Eddington limit, the higher (Lfir) / (L1450) found with these object (see Table 4) may 
suggest a higher ratio between SFR and SMBH accretion rate in these systems compared 
to the MAMBO detections among the optically bright z~6 quasars. This may indicate that 
the stellar bulge grows faster than the central SMBH in these faint/less massive quasar-host 
systems at z~6. 



6. Conclusions 

We present observations of millimeter and radio continuum and CO (6-5) line emission 
from the host galaxies of quasars at z~6. The new observations complete our MAMBO 
dust continuum survey of all the UV/optically faint z~6 quasars (mi45o > 20.2) discov- 
ered from SDSS. Combining with previous data, we calculate the average FIR and radio 
emission for a sample of 18 z~6 quasars with mi45o > 20.2. The mean FIR-to-AGN UV 
luminosity ratio of this faint quasar sample is about two times higher than that of the bright 
z~6 quasars (mi45o < 20.2). A fit to the average FIR and AGN bolometric luminosities 
of both the UV/optically faint and bright z~6 quasars, and the average luminosities of 
samples of submillimeter/millimeter-observed quasars at z~2 to 5, yields a relationship of 
Lfir ~ Lboi°'^^- The millimeter observations have detected the most FIR luminous objects 
among these faint z~6 quasars. The 250 GHz dust continuum detections in five of the 18 
faint quasars and CO (6-5) line detections in three of them (see Table 1 and 2) reveal dust 
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masses of a few 10^ M0 and molecular gas masses of 10^° M© in the quasar host galaxies, 
which are comparable to the dust and gas masses found in the bright z~6 quasars (Beelen 
et al. 2006; Bertoldi et al. 2003b; Carilli et al. 2007; Wang et al. 2008a, 2010). Their 
FIR luminosities are estimated to be 2.6 ~ 5.5 x 10^^ Lq, and the FIR and AGN bolo- 
metric luminosities follow a shallower relationship of Lpm ~ Lboi°'^^, in agreement with 
the starburst-AGN systems at lower redshifts. All these results suggest active star forma- 
tion in the host galaxies of the millimeter-detected and UV/optically faint quasars at z~6 
with SFRs of a few hundred MQyr~^. Moreover, the higher average FIR-to AGN UV lu- 
minosity ratios (Table 4) found with these objects and the shallow luminosity relationship 
suggest higher SFR-to-AGN accretion rate ratios than that of the more luminous/massive 
z~6 quasars. Further high-resolution imaging of the dust and molecular gas components 
(e.g., with ALMA) in the quasar host galaxies will be critical to constrain the host galaxy 
dynamical mass, black hole-spheroidal host mass ratio, gas and star formation rate surface 
density of these objects. 

This work is based on observations carried out with the Max Planck Millimeter Bolome- 
ter Array (MAMBO) on the IRAM 30m telescope, the Plateau de Bure Interferometer, and 
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many) and IGN (Spain). The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement by Associated Uni- 
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Humboldt Foundation through the Max-Planck-Forschungspreis 2005. Dominik A. Riechers 
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Fig. 1. — Left - VLA 1.4 GHz continuum image of the z=6.132 quasar J1319+0950. The 
beam size is 1.61" x 1.45",which is plotted on the bottom left, and the contour levels are (-3, 
-2, 2, 3, 4)xl5 /iJybeam"^. Right - PdBI 97 GHz continuum image of the quasar averaged 
over the line-free channels. The synthesized beam size (FWHM) is 3.5", and the contour 
levels are (-2, -1, 1, 2, 3) x 0.08 mJybeam"^. The cross in both plots denotes the optical 
position of the quasar (Mortlock et al. 2008). 



This preprint was prepared with the AAS E^TJtjX macros v5.2. 
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Table 1: Summary of the new observations 



Name 
(1) 



rcdshift 


mi450 


S250GHZ 


Sl.4GHz 


IAvco(6-5) 


FWHM 


Ocon 






mjy 


/ijy 


Jy kms^^ 


kms^^ 


mJy 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


5.7794±0.0008 


22.28 


2.37±0.49 


- 


0.37±0.07 


283±87 


0.14±0.04 


5.72 


20.97 


1.85±0.46 


195±22'' 


<0.32 


- 


<0.12 


5.82 


22.28 


1.29±0.54 


- 


- 


- 


- 


6.1321±0.0012 


19.65 


4.20±0.65 


64±17 


0.43±0.09 


537±123 


0.31±0.08 


6.12 


19.82 


l.OOiO.46'' 


23±18 


- 


- 


- 


5.92 


21.20 


0.09±0.63 


- 


- 


- 


- 


5.81 


21.65 


-0.36±0.61 


-28±18 


- 


- 


- 


5.87 


21.73 


-0.56±0.53 


-21±17 


- 


- 


- 


6.00 


21.77 


0.21±0.49 


- 


- 


- 


- 



SDSS J012958.51-003539.7 
SDSS J020332.39+001229.3 
SDSS J023930.24-004505.4 
ULAS J131911.29+095051.4 
CFHQS J150941. 78-174926.8 
SDSS J205321. 77+004706.8 
SDSS J214755.40+010755.0 
SDSS J230735.40+003149.0 
SDSS J235651. 58+002333.3 



Note - Column (1), name; Column (2), redshift. CO redshifts are presented for the two sources detected by the PdBI, and we 
list the redshifts measured with the UV quasar emission lines from the discovery papers for the other six sources (Jiang et al. 
2008, 2009); Column (3), rest-frame 1450 A magnitudes from the discovery papers (Mortlock et al. 2008, Jiang et al. 2008, 
2009). Column (4), MAMBO measurements of the 250 GHz dust continuum; Column (5), VLA measurements of the 1.4 GHz 
radio continuum; Column (6) and (7), line flux and FWHM of the CO (6-5) line emission; Column (8), measurement of the 
dust continuum at the CO (6-5) line frequency. The 250 GHz and 1.4 GHz detections arc marked as boldface. 
"Wang ct al. (2008b); ''Willott et al. 2007. 



Table 2: Previous observations of the 7721450 ^ 20.2 z~6 quasars 

Name 

(1) 

SDSS J000552.34-000655.8 
CFHQS J003311.40-012524.9 
SDSS J030331.40-001912.9 
SDSS J035349.76+010405.4 
NDWFS J142516.30+325409.0 
FIRST J142738.59+331242.0 
SDSS J163033.90+401209.6 
CFHQS J164121. 64+375520.5 
SDSS J205406.42-000514.8 
SDSS J231546.36-002357.5 
CFHQS J232908.28-030158.8 

Note -Reference: (1) Wang ct al. 2007; (2) Wang et al. 2008b; (3) Willott et al. 2007; (4) Wang et al. 2010; (5) McGreer et al. 
2006; (6) Bcrtoldi et al. 2003a. 

"The values are calculated from the absolute 1450 A magnitudes in the discovery papers (Cool ct al. 2006; McGreer et al. 
2006; Willott et al. 2007). The 250 GHz and 1.4 GHz detections are marked as boldface. 'The upper limits of LpiR are 



;dshift 


mi450 


S250GHZ 


Sl.4GHz 


LpiR 


Mdust 


Mgas 


Referei 






mJy 


A^Jy 


IOI2 L0 


10* M0 


IOIOM0 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


5.85 


20.83 


0.36+0.48 


40+130 


<3.4 


- 


- 


1,2 


6.13 


21.78" 


1.13+0.36 


-27+19 


2.6+0.8 


1.5 


- 


2,3 


6.08 


21.28 


0.23+0.51 


-85+62 


<3.5 


- 


- 


2 


6.05 


20.22 


1.20+0.46 


17+19 


<3.2 


- 


- 


2 


5.89 


20.62" 


2.27+0.51 


20+20 


5.4+1.2 


3.0 


2.0 


2,4 


6.12 


20.33" 


0.39+0.66 


1730+131 


<4.6 


- 


- 


2,5 


6.05 


20.64 


0.80+0.60 


14+15 


<4.2 


- 


- 


2,6 


6.05 


21.30" 


0.08+0.46 


-30+32 


<3.2 


- 


- 


2,3 


6.04 


20.60 


2.38+0.53 


17+23 


5.5+1.2 


3.1 


1.2 


2,4 


6.12 


21.34 


0.28+0.60 


31+16 


<4.1 


- 


- 


2 


6.42 


21.65" 


0.01+0.50 


14+22 


<3.4 


- 


- 


2 



derived with the 3cr upper limits of the 250 GHz flux densities. 
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Table 3: Derived parameters for the new observations 





Name 


LpiR 

1012 L© 


Mdust 

10* Mq 


Ll.4GHz 

LqHz-i 


^CO(6-5) 
10lOKkms-lpc2 


'^CO(l-O) 
10iOKkms-ipc2 


Mgas 

IOIOM© 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


J0129-0035 


5.2±0.9 


2.9 


- 


1.2±0.2 


1.5±0.3 


1.2 


J0203+0012 


4.4±1.1 


2.5 


0.118±0.013 


<1.0 


<1.3 


<1.0 


J0239-0045 


<3.9 


- 


- 


- 


- 


- 


J1319+0950 


10.0±1.3 


5.7 


0.045±0.012 


1.5±0.3 


1.9±0.4 


1.5 


J1509-1749 


<3.2 


- 


<0.038 


- 


- 


- 


J2053+0047 


<4.4 


- 


- 


- 


- 


- 


J2147+0107 


<4.4 


- 


<0.034 


- 


- 


- 


J2307+0031 


<3.8 


- 


<0.032 


- 


- 


- 


J2356+0023 


<2.8 


- 


- 


- 


- 


- 



Table 4: Average FIR and radio emission of the z~6 quasars 





Group 




Number" 


(Ll450) 

1012 L0 


(f250GHz) 
mJy 


< LpiR > 
1012 Lq 


(Lftr) 

(Ll450) 


Number"^ 


(fl.4GHz) 

MJy 


< Ll.4GHz >"* 

LqHz-i 


q 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 








Optically faint z~6 quasars 


with mi45o 


>20.2 








All objects 




18 


4.3 


0.85±0.12 


2.0±0.3 


0.46 


14 


19±6 


0.012±0.003 


1.63±0.15 


(radio quiet) 




16 


4.0 


0.79±0.13 


1.9±0.3 


0.46 


12 


2±6 


<0.013 


>1.57 


250 GHz detections 


5 


4.3 


1.85±0.20 


4.3±0.5 


1.01 


4 


45±10 


0.030±0.07 


1.59±0.11 


(radio quiet) 




4 


4.2 


1.86±0.23 


4.3±0.5 


1.03 


3 


1±12 


<0.024 


>1.68 


250 GHz non- 


■detections 


13 


4.3 


0.33±0.15 


<1.8 


<0.42 


10 


7±7 


<0.019 


- 


(radio quiet) 




12 


3.9 


0.32±0.15 


<1.8 


<0.46 


9 


2±7 


<0.015 


- 








Optically bright z 


~6 quasars 


with mi45o 


< 20.2 








All objects 




22 


17.6 


1.64±0.13 


3.8±0.3 


0.22 


22 


50±3 


0.033±0.002 


1.49±0.04 


(radio quiet) 




21 


16.9 


1.68±0.13 


3.9±0.3 


0.23 


21 


38±3 


0.026±0.002 


1.61±0.05 


250 GHz detections 


8 


19.1 


2.71±0.18 


6.3±0.4 


0.33 


8 


43±4 


0.029±0.003 


1.77±0.05 


250 GHz non- 


■detections 


14 


16.8 


0.55±0.18 


1.3±0.4 


0.08 


14 


59±5 


0.040±0.003 


0.94±0.15 


(radio quiet) 




13 


15.5 


0.58±0.18 


1.4±0.4 


0.09 


13 


32±5 


0.021±0.003 


1.23±0.15 



Note-We calculate the average flux densities using (f^) = J^ Wifi,^i/J^ Wi, Wi = l/fj;;, and {uv) = (X]™i) '^''' (Omont ct al. 
2003), where f^ ; and cr^i are the measurement and Icr rms for each object at 250 GHz or 1.4 GHz. For results with (f^) less 
than 3((T^), we take the stacking average value plus 3(0";^) as the upper limits of the average flux densities and calculate the 
corresponding upper limits for the luminosities. 

" Number of sources observed at 250 GHz. ^ Average luminosity at rest-frame 1450 A, derived from mi45o. '^ Number of 
sources observed at 1.4 GHz. '^ A radio spectral index of -0.75 (Condon 1992) is adopted here to calculate the rest frame 1.4 
GHz lunimosity. 
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Fig. 2. — The CO line spectra (left) and velocity-integrated images (right) of the three 
z~6 quasars observed with the PdBI. The continuum emission underlying the CO spectra 
is removed for the two detections, J1319+0950 and J0129— 0035. The synthesized beam 
size (FWHM) of the CO velocity-integrated image of J1319+0950 is ~3.5", and ~ 5" for 
J0129— 0035 and J0203-I-0012. The top abscissae of the CO spectra give the redshift range 
of the spectral windows, and zero velocity in each panel corresponds to the optical/near- 
infrared redshifts of z=6.127, 5.780, 5.720 from the discovery papers (Mortlock et al. 2008; 
Jiang et al. 2009). The solid lines in the spectra show Gaussian fits to the line emission. 
The CO peak positions are consistent with the optical quasars for the two detections. 
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Fig. 3. — The FIR versus radio luminosity plot of the MAMBO-detected radio-quiet z~6 
quasars that have radio 1.4 GHz continuum observations. The filled blue and red stars 
are the average values for the MAMBO-detected optically bright (mi45o < 20.2) and faint 
(nii45o > 20.2) z~6 quasars, respectively. The open stars represent the stacking average or 
upper limits of the non-detections in the two samples. Only the radio-quiet objects in the 
z~6 quasar samples are plotted here. The crosses are a sample of IRAS star forming galaxies 
from Yun et al. (2001). The dashed line shows the average radio-FIR correlation in star 
forming galaxies (q = log (Lfir/3.75 x 10^^ L0) — log (Li.4ghz/L0 ■ Hz~^) = 2.34, Yun et al. 
2001). The dotted lines represents q values five times higher and lower, which inlcudes 98% 
of the IRAS galaxies. 
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Fig. 4. — The FIR and AGN bolometric luminosity correlations of the z~6 quasars, together 
with samples of local optically selected IR luminous quasars hosted by ULIRGs (Hao et al. 
2005; Zheng et al. 2002) and the average luminosities of the (sub)mm-observed quasars at 
redshift z~2 to 5 (Omont et al. 2001, 2003; Priddey et al. 2003a; Carilh et al. 2001). The 
filled squares and circles are the MAMBO detections in the bright (mi45o < 20.2) and faint 
(mi45o > 20.2) z~6 quasar samples, respectively, and the arrows indicate 3cr upper limits 
of the non-detections (see Section 3.2). The z=6.0 quasar, SDSS J130608. 26+035626. 3, was 
not detected at 250 GHz (Wang et al. 2007), but was detected by SCUBA at 350 GHz 
(Priddey et al. 2003b). We estimate the FIR luminosity of this object with the SCUBA 
flux density and plot it as filled square. The filled triangles and diamonds are the mean FIR 
luminosities for the (sub)mm detections in each sample, derived by averaging flux densities 
at 250 GHz or 350 GHz. The open symbols represent the mean luminosities averaged with 
both the millimeter detections and non-detections in each sample. The dotted line is a 
fit to the average luminosities (open symbols) of all the high-z samples, i.e. log(LFiR) = 
0.621og(Lboi) + 3.9, and the dashed line shows the fit to the submillimeter or millimeter 
detected sources in all high-z samples and the local ULIRGs, log(LFiR) = 0.451og(Lboi) + 6.6. 



